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Abstract 

The equilibrium thermodynamic properties of the SU(iV) plasma at finite temperature 
are studied non-pcrturbatively in the large-iV limit, via lattice simulations. We present high- 
precision numerical results for the pressure, trace of the energy-momentum tensor, energy 
density and entropy density of SU(iV) Yang-Mills theories with N = 3, 4, 5, 6 and 8 colors, 
in a temperature range from 0.8T C to 3.4T C (where T c denotes the critical dcconfincmcnt 
temperature). The results, normalized according to the number of gluons, show a very mild 
dependence on N, supporting the idea that the dynamics of the strongly-interacting QCD 
plasma could admit a description based on large- N models. We compare our numerical data 
with general expectations about the thermal behavior of the deconfined gluon plasma and 
with various theoretical descriptions, including, in particular, the improved holographic QCD 
model recently proposed by Kiritsis and collaborators. We also comment on the relevance of an 
AdS/CFT description for the QCD plasma in a phenomenologically interesting temperature 
range where the system, while still strongly-coupled, approaches a 'quasi-conformal' regime 
characterized by approximate scale invariance. Finally, we perform an extrapolation of our 
results to the N — > oo limit. 



PACS numbers: 12.38.Gc, 12.38.Mh, 11.15.Ha, 25.75.Nq, 11.10. Wx, 11.15.Pg, 11.25.Tq 



1 Introduction 



A direct implication of asymptotic freedom in non-Abelian gauge theories [1,2] is that hadronic 
matter should undergo a change of state to a deconfined phase at sufficiently high temperatures 
or densities [3,4]. This theoretical prediction has been subject to extensive experimental investi- 
gation through heavy- ion collisions for more than two decades. The evidence obtained from the 
SPS and RHIC experiments led to the conclusion that 'a new state of matter' was created [5-9], 
which reaches rapid thermalization, and is characterized by very low viscosity values, making it 
a nearly ideal fluid [10]; a brief summary of experimental results can be found in ref. [11]. 

On the theoretical side, however, the understanding of such strongly-interacting quark-gluon 
plasma (sQGP) is still an open issue [12]. Relativistic fluidodynamics — see ref. [13] for an intro- 
duction to the subject — provides a successful phenomenological description of the experimental 
data, but it is not derived from the first principles of QCD. On the other hand, the weak-coupling 
expansion in finite-temperature QCD is non-trivial: it contains terms proportional to odd powers 
of the gauge coupling [14], as well as contributions from diagrams involving arbitrarily large num- 
bers of loops [15, 16], and is slowly convergent; finally, it shows quite large discrepancies with the 
results determined non-perturbatively from lattice simulations at temperatures close to the de- 
confinement [17-22]. One possibility is that this mismatch is due to intrinsically non-perturbative 
contributions proportional to T 2 [23,24], which could perhaps be modelled in terms of a modified 
bag model [25] with finite thickness; related ideas have also been discussed in refs. [26-30]. Alter- 
native descriptions of the sQGP phenomenology are based on various types of effective degrees 
of freedom [12,31-35]. 

For a first-principle description of the QCD plasma, however, the numerical approach on 
the lattice is essentially the main tool. In particular, the lattice determination of equilibrium 
thermodynamic observables in the SU(3) gauge theory is considered as a solved problem [36]; on 
the other hand, several studies of finite-temperature QCD with dynamical fermions have begun 
appearing during the last decade [37-44], and today large-scale unquenched simulations can be 
done at or close to the physical point; recent results are reviewed in refs. [45-47]. 

Yet, studies of the thermodynamics of the pure-gauge sector can still be relevant from a 
fundamental perspective, as they capture the essential qualitative features of the deconfinement 
phenomenon, and offer several advantages. 

Firstly, the deconfining transition at finite temperature in pure Yang-Mills (YM) models is 
characterized by a well-defined theoretical setup, with a clean symmetry pattern: in the Euclidean 
formulation, the action is invariant under center symmetry, and the trace of the Wilson line (or 
Polyakov loop) along the compactified Euclidean time direction serves as the order parameter 
distinguishing the confined (symmetric) from the deconfined phase (with spontaneously broken 
symmetry). By contrast, in full QCD with physical quark masses the center symmetry is explicitly 
broken by the quarks, and no unambiguous definition of a critical temperature exists, with the 
change of state from the confined to the deconfined phase being a cross-over. 

Secondly, calculations in lattice QCD including fermionic degrees of freedom are computation- 
ally much more demanding than in YM models, and all presently known methods to implement 
the fermionic sector on the lattice have advantages and drawbacks [48]. On the other hand, 
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numerical simulations in the pure gauge sector can attain much higher precision, and allow to 
detect fine details in the behavior of various observables. 

Finally, the study of SU(iV) gauge models with A > 3 colors is relevant for the large-A 
limit. Historically, the idea of studying QCD in the large-A limit was first put forward in 1974 
by 't Hooft [49], who proposed to consider 1/N as an expansion parameter: the limit for N 
going to infinity at fixed bare 't Hooft coupling Ao = g^N (go being the bare YM coupling) 
yields a simpler theory, characterized by dominance of a subset of Feynman diagrams (which, 
in a double-line notation, can be represented as planar ones), and amenable to a topological 
expansion in terms of surfaces of increasing genus. Indeed, the simplifications taking place in the 
large- A model can be useful to obtain insight into certain non-trivial non-perturbative features of 
QCD [50,51]. As finite- N corrections are of order 0(N~ 1 ) for QCD (or 0(A~ 2 ) for pure YM), one 
may expect that predictions derived from large- N QCD have phenomenological relevance; several 
issues related to this expectation have been investigated on the lattice by different groups — see 
refs. [52-54] and references therein. 

However, from a more modern perspective, the large- N limit has much farther-reaching im- 
plications (see, e.g., refs. [55,56] and references therein), and also entails deep connections with 
string theory. Indeed the large- N limit is one of the crucial ingredients underlying the AdS / CFT 
correspondence conjectured by Maldacena [57], according to which the large- A limit of the max- 
imally supersymmetric N = 4 supersymmetric YM (SYM) theory in four dimensions is dual to 
type IIB string theory in a AdS§ x S" 5 space. In principle, this highly non-trivial correspondence 
allows one to study the strongly-coupled regime of field theory by mapping it to the weak-coupling 
limit of a gravity model. The AdS/CFT correspondence is an example of a duality relating the 
large- A limit of gauge theories with U(A) structure group and string theory. The mathematical 
meaning of such dualities can be summarized in the following terms: gauge theories deal with 
invariants originating from moduli spaces of connections on principal bundles, whereas string 
theory studies the invariants arising from spaces of maps from certain classes of domains onto 
different targets; large- A dualities relate the generating functions of such invariants, i.e. the 
partition functions of these theories [58] . 

Although the Maldacena conjecture was originally formulated for the N = 4 SYM theory, 
the techniques underlying the AdS/CFT correspondence are often applied to theories which 
are expected to have properties more similar to QCD: in the so-called 'top-down' approach, 
this is achieved by explicitly breaking some symmetries of the N = 4 theory. This allows to 
build models with a non-trivial hadron sector, which can be studied non-perturbatively using 
the dual gravity formulation [59]; for the meson sector, the predictions obtained from these 
calculations have been compared with (quenched) spectroscopy calculations in SU(A) lattice 
gauge models [60,61], finding satisfactory agreement. Another field where analytical techniques 
inspired by the AdS/CFT correspondence are applied is in the description of hydrodynamic and 
thermodynamic properties relevant for a strongly interacting system, like the QCD plasma — see 
refs. [62-70] and references therein. In principle, the appropriate parameters for the gravity 
models dual to the sQGP can be pinned down using the combined constraints coming from 
experimental data and lattice QCD simulations [71]. 

On the other hand, similar analytical techniques have also been used in a 'bottom-up' ap- 
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proach, building five-dimensional holographic models that should reproduce the non-perturbative 
features of QCD [72-75]; this approach has proven quite successful for the hadron properties, but 
is not completely satisfactory for the QCD plasma thermodynamics [76]. More recently, a related 
approach has been used to build an "improved holographic QCD model" [77-80] through a dual 
description in terms of a dilaton-gravity system in five dimensions; this model reproduces the 
main non-perturbative features of the SU(3) YM spectrum as determined in lattice calculations, 
as well as the equilibrium thermodynamic properties of the sQGP [80-82]; furthermore, the model 
can also be used to derive predictions for observables which are more difficult to study on the 
lattice, such as the plasma bulk viscosity, drag force and jet quenching parameter [83]. Related 
works include refs. [84-93]. 

Clearly, all models aiming at a description of N = 3 QCD (or YM) in terms of predictions 
based on the large-iV limit implicitly rely on the assumption that the features of the N = 3 
theory are 'close enough' to those of its N = oo counterpart. While a priori this assumption is 
not guaranteed to be true, there is strong numerical evidence that this is indeed the case [52] . In 
this context, the aim of the work reported here was to study non-perturbatively the equilibrium 
thermodynamics properties at finite temperature in SU(iV) YM theories with N = 3, 4, 5, 6 and 
8 colors, via lattice simulations. Related studies include works by Bringoltz and Teper [94, 95] 
and by Datta and Gupta [96-98]. 

The structure of this paper is the following: in sec. El we briefly review some theoretical 
expectations for the sQGP thermodynamics, which can be compared with lattice results at large 
N. The setup of our lattice calculation is formulated in sec. [31 while in sec. [5] we present the 
numerical results for the pressure, trace of the stress-energy tensor, energy and entropy densities, 
and compare them with the theoretical expectations; we discuss the main uncertainties involved 
in our calculation, and the extrapolations to the continuum and to the large- N limit. Finally, we 
conclude in sec. [5] with a summary of our findings, a discussion about their implications, and an 
outlook on future work. Preliminary results of this study have been presented in ref. [99] . 

2 Pictures of the Unicorn: Theoretical descriptions of the sQGP 
thermodynamics 

In the literature, there exist many theoretical approaches to describe the sQGP thermodynamics; 
the elusiveness of the QCD plasma, which has inspired the comparison with the mythological 
creature (see ref. [100] and ref. [3] therein), is not only a matter of experimental issues, but is 
also related to the theoretical difficulties in working out a quantitatively successful description 
for a strongly-interacting system. 

While a review of the theoretical approaches would be a task beyond the scope of this paper 
(we refer the reader to reviews available in the literature, like, e.g., refs. [12,31]), in this section 
we just recall certain theoretical expectations, for which a comparison with the numerical results 
of N > 3 simulations may be of interest. 

Roughly speaking, the logical organization of this section is from the perturbative to the 
non-perturbative, from QCD to strings, from lower to higher dimensions. More precisely, we 
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start by briefly discussing the key issues in the weak-coupling expansion for thermal QCD, and 
the possibility to resolve its mismatch with lattice results at temperatures of the order of the 
deconfinement temperature T c by including terms of non-perturbative origin, which could perhaps 
be described in terms of a "fuzzy bag model". Then we move to the description of the sQGP 
by gravity models, reviewing the basic features of the improved holographic QCD model, and 
recalling some AdS/CFT predictions which can be compared with our numerical results. 

2.1 Weak-coupling expansions in hot QCD 

Due to asymptotic freedom, the weak-coupling expansion is expected to be a natural analytical 
approach to the thermodynamics of hot QCD [32]. In particular, at very high temperatures 
QCD develops a hierarchy of momentum scales, with a clear separation between 'hard', 'soft' and 
'ultrasoft' modes, whose parametric dependence on the temperature T and on the coupling g is 
0(T), O(gT) and 0(g 2 T), respectively. By integrating out the hard modes, one can derive a 
dimensionally-reduced effective theory [101-104], whereby the QCD pressure can be studied in a 
weak-coupling expansion [20,105,106] and some of the intrinsically non-perturbative terms [15,16] 
can be estimated numerically [22, 107-109]. However, the convergence of the perturbative series 
appears to be rather slow, and the most recent results (for the pure gauge sector) [22] exhibit large 
deviations from the corresponding lattice calculations at temperatures 0{T C ) [36]; perhaps the 
situation could be improved using a center-symmetric dimensionally-reduced theory [110-112]. 

2.2 Non-ideal contributions to the QCD pressure and a "fuzzy bag" model 

The mismatch between the weak-coupling expansion for the pressure in thermal QCD and lattice 
results is particularly severe for T c < T < 3T C , see, e.g., fig. 6 in ref. [22]. Pisarski observed that in 
this temperature range the lattice results for the trace of the stress-energy tensor (often denoted 
as A) in the SU(3) gauge theory appear to be well-described by a T 2 behavior [23,24], which 
is not captured by the weak-coupling expansion. For this reason, he argued that the expression 
of the pressure p, as a function of the temperature in the range between the maximum of A/T 4 
(which, for SU(3), is located around 1.1T C ) and the onset of the perturbative behavior, should be 
of the form [24] : 

P{T) = fpcrtT 4 - B [uzzy T 2 - B MlT + . . . (1) 

where / pert includes the perturbative contributions, -B fuzzy is a "fuzzy bag" term, and B MIT is the 
term associated with the usual MIT bag model [25]. The idea behind the B fuzzy term is that 
of a QCD vacuum where deconfined "bubbles" have a surface of finite thickness; if the degrees 
of freedom associated with fluctuations of this thick bag surface are collective ones and involve 
the whole surface width, then they are expected to be difficult to excite, thus they should not 
lead to (unphysical) multiplicities in the zero-temperature hadron spectrum. On the other hand, 
assuming that the surface of the fuzzy bag extends over a width between 0.3 fm and the typical 
hadron size of 1 fm would make it plausible that it should become relevant for temperatures in a 
range between T c and 3T C . 
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From this point of view, one can consider the fuzzy bag model as a phenomenological descrip- 
tion for the QCD vacuum and for the plasma thermodynamics; such description involves two free 
parameters (i? fuzzy and B MIT ). One possible way to relate this description to more fundamental 
degrees of freedom is in terms of an effective theory of Wilson lines, whose dynamics could be 
described by some random matrix model [113,114]; in this respect, confinement would arise from 
the eigenvalue repulsion in the large- iV limit — an issue which has been studied in a finite volume 
in ref. [115]. 



2.3 Thermodynamics from the improved holographic QCD model 



A radically different description for the YM spectrum and for the deconfined gluon plasma can 
be obtained from the improved holographic QCD (IHQCD) model proposed by Kiritsis and 
collaborators [77-83], which is an AdS/QCD model based on an Einstein-dilaton gravity theory 
in five dimensions, where the fifth (radial) direction is dual to the energy scale of the SU(iV) 
gauge theory. On the gravity side, the bulk fields are the metric, the dilaton and the axion, which 
are respectively dual to the energy- momentum tensor, to the trace of F 2 , and to the trace of FF 
on the gauge side. The action is defined as: 



Sihqcd 



-M P N' Z I dPxy/g 



R 



^f + V{\) 



+ 2MpN 2 [ d A xVh K, 
JdM 



(2) 



where <fr is the dilaton field, while A = exp($) is identified with the running 't Hooft coupling of 
the dual SU(iV) YM theory, Mp denotes the five-dimensional Planck scale, and K is the extrinsic 
curvature of the boundary a la Gibbons-Hawking. Note that the effective five-dimensional Newton 
constant is G5 = 1/ (lQitMpN 2 ^ , which becomes small in the large- N limit. The dilaton potential 
V(X) is defined according to the requirements of asymptotic freedom with a logarithmically 
running coupling and linear confinement in the ultraviolet (UV) and infrared (IR) limits of the 
dual gauge model; one possible Ansatz is: 



V(X) 



12 



1 + V0X + V1A 4 / 3 Jbg (1 + y 2 A 4 /s + V3A 2 ) 



(3) 



where I is the AdS scale (which only defines the overall normalization of the potential). Vq, V\, V2 
and V3 are free parameters: two of them can be fixed by imposing that the dual model reproduces 
the first two (scheme-independent) perturbative coefficients of the SU(iV) /^-function, and one is 
left with two independent parameters. The latter can be adjusted by requiring that the values 
for two independent physical quantities in the dual model match the corresponding YM results 
obtained from lattice simulations. 

In this model, the first-order transition from a confined to a deconfined phase (respectively 
dual to a thermal-graviton- and to a black-hole-dominated regime of the five-dimensional gravity 
theory) is associated with the appearance of a non-trivial gluon condensate. 

Strictly speaking, the IHQCD model is expected to hold in the large-iV only, while at finite 
N one expects corrections; in particular, the calculations in the gravity model neglect string 
interactions, which are expected to become important above a cut-off scale MpiV 2 / 3 ; for the 
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values of the parameters used in ref. [82], this would approximately correspond to 2.5 GeV in 
SU(3)— see also ref. [116]. 

2.4 A prediction from J\f = 4 SYM thermodynamics: The entropy deficit 

While the IHQCD model discussed in subsec. 12.31 provides a quantitatively accurate descrip- 
tion — rather than a prediction — for the non-perturbative features of the strongly-interacting YM 
spectrum and of the deconfined plasma, and is based on a non-supersymmetric, non-critical string 
formulation [117,118], a more fundamental/less phenomenological approach consists in deriving 
general properties of the sQGP directly from the AdS/CFT setup — see refs. [66-70] and references 
therein. 

It is worth stressing that, in principle, there are no obvious reasons to expect the large-A 
limit of M = 4 SYM to be a theory "close" to real-world QCD: as opposed to the latter, the 
former is maximally supersymmetric, conformally invariant, lacks confinement and spontaneous 
breaking of chiral symmetry. The field content of the two theories is also very different, given 
that SYM features fermions in the adjoint, rather than in the fundamental, representation, and 
that the number of colors N is infinite, rather than 3 as in QCD. 

Yet, this approach has led to the analytical derivation of important results, with potential 
relevance also for phenomenology. A celebrated example is the value for the shear viscosity to 
entropy density ratio rj/s = l/(4-7r) [62], which is expected to hold for all gauge theories in the 
limit of infinite 't Hooft coupling [119], and which is conjectured to be a universal lower bound 
for all substances existing in Nature [120]. 

In the context of the equilibrium thermodynamics quantities that we studied in the present 
work, one important prediction from large- N J\f = 4 SYM concerns the entropy density s, which 
can be evaluated in the limit of infinite coupling A — > oo by identifying it with the Bekenstein- 
Hawking entropy density of the dual gravity theory: the result is found to be equal to 3/4 of the 
entropy density calculated in the limit of zero coupling (denoted as so) [121]. This result was 
later refined in ref. [122], with the evaluation of the leading correction at large but finite A values: 

- = z + |c(3)(2Ar 3 / 2 + ... 

so 4 32 

where £ is the Riemann function, and the dots denote subleading corrections; as Apery's constant 
C(3) ~ 1.2020569 is a positive number, the asymptotic value of s/sq is approached from above 
when A tends to infinity. 

Note that — as opposed to QCD or to the IHQCD model discussed in subsec. 12.31 — in the 
supersymmetric model A is merely a parameter of the theory, which does not run with the 
energy. 

3 Lattice formulation 

In this work, we performed a non-perturbative study of SU(iV) YM theories with A = 3, 4, 5, 6 
and 8 colors, regularizing them on a four-dimensional Euclidean hypercubic, isotropic lattice with 
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periodic boundary conditions in all directions, and estimating the expectation values of physical 
observables via Monte Carlo simulations. In the following, a denotes the lattice spacing, N s (Nt) 
is the number of lattice points along each of the three space-like (along the time-like) directions, 
while V = (aN s ) 3 and T = \/{aNt) denote the spatial volume and temperature, respectively. 
The degrees of freedom of the theory are represented as U^(x) G SU(iV) matrices of size N x N, 
defined on the oriented elementary bonds (x, p) of the lattice; their dynamics is governed by the 
standard (unimproved) Wilson gauge action Sw [123]: 



where the summation ranges over all the elementary lattice plaquettes p, with U p denoting the 
path-ordered product of the link variables around each plaquette. The partition function Z = 



where DU denotes integration over the Haar measure for all U^(x), and (3 is related to the bare 
lattice gauge coupling go through: = 2N/gQ. The expectation values of the various observables 
are estimated via a dynamic Monte Carlo simulation sampling the space of configurations of 
the system. The configuration update algorithm to generate the Markov chain is based on a 
combination of heat-bath steps for the SU(2) subgroups of SU(iV) [124-126] and full SU(iV) 
overrelaxations [127,128] — see also ref. [129] — , except for a subset of the T = simulations, which 
were performed using the Chroma library [130] on the QCDOC machine [131] in Regensburg. 

The physical scale for the SU(3) gauge group was set using the values of Sommer's parameter 
ro [132] taken from ref. [133]. For the other SU(iV > 3) groups, the physical scale was defined 
from interpolation of high-precision measurements of the string tension a [134, 135]; for the (5 
values beyond the range where high-precision string tension measurements are available, this 
was combined with a three-loop perturbative expansion [136] in terms of the Parisi mean-field 
improved bare lattice coupling [137], taking lattice corrections into account, following ref. [138]. 

The parameters of the main set of lattice simulations performed in this work (denoted as 
'ensemble A' in the following) are summarized in tab. Q] 

We focus on the following equilibrium thermodynamic quantities: free energy density (/), 
pressure (p), energy density (e), trace of the energy-momentum tensor (A) and entropy density 




p 



(4) 



Z(T, V) is given by: 




(5) 
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N 


N s 


N t 


/^min 




Np 


A/3 


/3-range 


T = statistics 


finite-T statistics 


3 


20 


5 


5.7 


6.6 


181 


0.005 


[5.7,5.795] 


3000 


12000 
















[5.8,5.845] 


5800 


33600 
















[5.85,5.895] 


5800 


21600 
















5.9 


6400 


21600 
















[5.905,6.2] 


5800 


21600 
















[6.205,6.6] 


3000 


12000 


4 


16 


5 


10.5 


11.8 


131 


0.01 


[10.5,11.5] 


16000 


24000 
















[11.51,11.67] 


4200 


12300 
















11.68 


7400 


12300 
















[11.69,11.8] 


4200 


12300 


5 


16 


5 


16.75 


18.45 


35 


0.05 


[16.75, 18.45] 


2500 


7500 


6 


16 


5 


24.4 


26.9 


26 


0.1 


[24.4, 25.5] 


10000 


20000 
















[25.6,26.9] 


5200 


8400 


8 


16 


5 


43.9 


48.1 


43 


0.1 


[43.9,48.1] 


3200 


9200 



Table 1: Parameters of the main set of lattice simulations used in this work (ensemble A). iV 
denotes the number of colors, N s and Nt are the space-like and, respectively, time-like sizes of 
the system in lattice units, Np is the number of /3-values that were simulated, separated by A/3 
from each other, in the /3 min < < /? max interval. For each subset of /3-values in this range, the 
T = and finite-T statistics at each (3 are shown in the last two columns. 



(*): 



/ 
P 
e 

A 

s 



T 

T logZ 
, dlogZ 



dV 
T 2 dlogZ 
~V 



dT 



v 



3p 



-\ogZ + -. 
V T 



(6) 
(7) 

(8) 
(9) 
(10) 



In the thermodynamic limit one has: p = —f, while in a finite hypertorus in general this relation 
gets infrared corrections [139-144]. In particular, the infrared corrections to the equation of state 
of a gas of free gluons in a finite volume were derived analytically in ref. [142]; such corrections 
could be particularly relevant for lattice simulations at very high temperatures [145], especially 
for the case of a fixed value of the aspect ratio of a time-like cross-section of the system LT. 
Although the investigation of such a limit (i.e. at a very fine lattice spacing a, with a shrinking 
physical space-like volume) is interesting in its own [146, 147] , in the present work we restricted 
our attention to relatively low temperatures, in a range up to about 3T C , where — due to screen- 
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ing [16] — such effects are not visible within the level of precision of our simulations. This was 
investigated in our preliminary study reported in ref. [99], which showed that, for the setup of 
the calculations presented in this work, finite-size effects are negligible. 

In the regime where finite-volume corrections can be neglected, the pressure — relative to its 
T = vacuum value — can be conveniently estimated by means of the "integral method" [148]: 

p ~ — losZ = —. — 5 — / dp — - , , (11) 

p V S a^NlNtJfo ^ d(3' ' { ' 

where (3q is chosen to be a point sufficiently deep in the confined phase (corresponding to a 
temperature at which the pressure differs from its T = value by a negligible amount). 
The dimensionless ratio p/T A can be readily evaluated from: 

J^ = 6j\(3'((U p ) T -(U p ) ), (12) 

where the notation (■ ■ - )t indicates the expectation value in the system at temperature T (which 
is estimated from simulations on a lattice of sizes N^ x N ), while (• • • )o denotes the expectation 
value at zero temperature (that can be estimated from measurements on a lattice of size iV 4 ). 

The integration in eq. (|12p was performed numerically, dividing the integration range in (5 
in N/3 — 1 equally spaced intervals (see tab. []]) and according to the methods discussed in the 
Appendix of ref. [149]. 

A/T 4 can be expressed as: 

— = 6N? ((U p )o - (U P ) T ) ■ a£. (13) 

The energy and entropy densities can then be obtained indirectly, as linear combinations of 
p/T and A/T 4 , namely: 

A + 3P (14) 



rp4 rp4 

and: 

s A + Ap 



(15) 



2^3 ^4 

Finally, it should be noted that lattice cut-off effects induce a systematic deviation from the 
continuum values of these quantities. In the infinite-temperature limit, such effect is quantified by 
a dimensionless ratio denoted as Rj(Nt), which can be evaluated by numerical integration [150]. 
In the present work, we always include this factor in the definition of the lattice Stefan-Boltzmann 
(SB) limit for the various thermodynamic observables. 



4 Results 

In this section, we present the numerical results from our main set of simulations on lattices with 
Nt = 5 (ensemble A) and compare them with the theoretical predictions in subsection 14.11 The 
systematic uncertainties of our calculations are discussed in subsection 14.21 where we also present 
an extrapolation of the thermodynamic observables to the large- limit. 
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Figure 1: (Color online) The dimensionless ratio p/T 4 , normalized to the lattice SB limit it 2 (N 2 — 
l)i?/(A r t)/45, versus T/T c , as obtained from simulations of SU(iV) lattice gauge theories on Nt = 5 
lattices. Errorbars denote statistical uncertainties only. The results corresponding to different 
gauge groups are denoted by different colors, according to the legend. The yellow solid line 
denotes the prediction from the improved holographic QCD model from ref. [82] (with a trivial, 
parameter-free rescaling to our normalization). 



4.1 Results from the simulation ensemble A 

The data obtained from our high-precision simulations reveal very strong qualitative and quantita- 
tive similarities between the results (normalized to their lattice SB limits) for the thermodynamic 
observables, in all of the gauge groups investigated in agreement with the findings of refs. [94,96]. 
This is manifest in figs. [H El E] and [H where we show the dimensionless ratios p/T 4 , A/T 4 , e/T 4 
and s/T 3 , normalized to their respective SB limitsU are plotted against T/T c . 

In these plots, we also show the comparison with the curves obtained from the improved 
holographic QCD model in ref. [82], assuming eq. ([3j) as an Ansatz for the dilaton potential, and 
by fitting two free parameters via a comparison with lattice results taken from other works. 

1 This feature has an interpretation in the quasi-particle model discussed in refs. [151,152]. 
2 Since the SB limit of A is zero, we conventionally normalize A/T 4 to the SB limit of p/T 4 . 
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Trace of the energy-momentum tensor 
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Figure 3: (Color online) Same as in fig. [TJ, but for the e/T 4 ratio, normalized to the SB limit. 



Note that in ref. [82], by looking at the comparison with the N = 3 lattice results for A/T 
from ref. [36], it was pointed out that the slight discrepancy between the improved holographic 
QCD model and the lattice results in the region of the peak (which for the SU(3) lattice data 
is located at T ~ 1.1T C , and is slightly lower than the IHQCD model curve) was likely to be a 
finite-volume lattice artifact. Here, fig. [2] indeed confirms this: our results for the SU(3) gauge 
group are consistent with ref. [36], while the A/T 4 maximum for the N > 3 gauge groups is 
higher and located closer to T c . This is related to the fact that the deconfining phase transition, 
which is a weakly first-order one for SU(3), becomes stronger when N is increased [153-155] — see 
also refs. [156-158] — and correspondingly the value of the correlation length at the critical point 
becomes shorter. 

In order to study the relevance of AdS/CFT effective models for the sQGP, it is interesting 
to investigate 'how close' to being conformal the deconfined system is. To address this issue, in 
fig. [5] we plot the lattice equation of state in a (e,p) plane (in a format similar to analogous plots 
in refs. [96, 159], up to a different normalization of the horizontal axis), where the conformal limit 
is described by the straight line through the origin and the point corresponding to the SB limit in 
the top-right corner. As the temperature is increased, the points tend to approach the conformal 
line from below. As it has been pointed out in ref. [96] (by looking at results of SU(3) and SU(4) 
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p{z) equation of state and approach to conformality 
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Figure 5: (Color online) The SU(iV) lattice equation of state, expressed as p{e), and the approach 
to conformality (dotted line). The dashed curve denotes the prediction from the SU(3) weak- 
coupling expansion [22]. 

gauge groups), this behavior favors perturbation theory [22] over a conformal description; our 
data, displayed in the same format as in ref. [96], appear to confirm this result. However, at 
temperatures about 3T C , although both e and p are still far from their SB values, the points (for 
all groups) lie close to the conformal line, so it is reasonable to expect that an AdS/CFT model 
could also provide a good effective description at these temperatures. 

In this respect, one observable which is particularly interesting to consider is the entropy 
density s, as in the literature it is often mentioned that the relative entropy deficit (with respect 
to the non-interacting limit) that is observed in lattice QCD simulations is close to the one 
predicted by AdS/CFT for the strongly-coupled limit of M = 4 SYM, which is 1/4 [121]. In 
fact, AdS/CFT also predicts the leading dependence of the s/so ratio (with so denoting the 
entropy density in the free case) on the 't Hooft coupling A of the supersymmetric theory, at 
the first perturbative order around the strong-coupling limit [122]. However, this information 
cannot be directly used for a comparison with the YM model, because the 't Hooft coupling of 
the supersymmetric model cannot be identified with its counterpart in the non-supersymmetric 
gauge model: the former does not run with the energy and is simply a parameter of the theory, 
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while the latter is a renormalized coupling which runs with the momentum scale and, in general, 
is scheme-dependent. One possible strategy to bypass this potential ambiguity and try to work 
out a comparison between the super symmetric model and the YM theory could be by considering 
the screening masses [160]. 

In general, insisting on identifying the A parameter of N = 4 SYM with some particular 
definition of the running 't Hooft coupling in the YM theory would be a scheme-dependent, 
and thus ambiguous, procedure. However, one may still argue that one could get at least some 
general, semi-quantitative information about how well the supersymmetric model captures the 
behavior of the YM theory (in the regime where it becomes 'quasi-conformal') by identifying the 
supersymmetric theory parameter with the YM running coupling evaluated in some physically- 
motivated, 'good' scheme (meaning a scheme for which low-order perturbation theory calculations 
work well down to energies close to the non-perturbative scale of the theory). 

In this context, we mention that in the recent literature there are various studies discussing the 
running of the coupling in large- N lattice gauge theories, with different approaches and in different 
schemes [138, 161-163], pointing to a mild, well-behaved dependence on N . In particular, one 
popular choice to define the running coupling is the Schrodinger functional (SF) scheme [164,165]; 
the behavior of the running coupling in the SF scheme in the SU(4) lattice gauge theory has been 
investigated by Lucini and Moraitis [163]: they found good agreement with two-loop perturbation 
theory down to energies comparable to the square root of the string tension. Assuming that such 
agreement also holds in the large- N limit [138, 163], one can estimate the running coupling for all 
gauge groups considered in the present work, and convert between different schemes — including, 
in particular, the modified minimal-subtraction (MS) scheme — using the expressions for Arg^/^fa 
and T c j \fa at different iV which can be taken from the literature [138, 155]. Following this 
approach, one can obtain a plot of the entropy versus an approximate estimate of the running 
't Hooft coupling in the MS scheme, which is shown in fig. [6j at temperatures around 3T C (where 
the plasma approaches a 'quasi-conformal' regime), the physical 't Hooft coupling for SU(3) is 
around 5.5 (in agreement with ref. [36]), and still large for all SU(iV) groups. 

Incidentally, it may be interesting to note that, if one were to directly identify the A parameter 
of the N = 4 model with the YM running 't Hooft coupling evaluated in the MS scheme, then 
eq. ([5]) would predict a curve (dashed line in fig. [6]) that is smoothly approached by the lattice 
results in the temperature regime where the AdS/CFT description is expected to begin to work 
(right end of the plot). Also, the values of the 't Hooft coupling in that regime agree with the 
estimate obtained in ref. [166] from the comparison of the drag force on heavy quarks in the 
N = 4 model and in lattice QCD simulations with dynamical quarks, and appear to be roughly 
comparable with the lower bound of the window for this parameter identified in ref. [71], namely: 
A ~ 10-25 (although the gravity model discussed there also includes a Gauss-Bonnet term [167], 
and the constraints on the parameters are derived by comparing with experimental data and with 
results of lattice QCD simulations with dynamical fermions [39-42]). 

However, we stress once more that the comparison suggested by fig. [6] should be taken cum 
grano salts. First and foremost, because of the ambiguity and scheme-dependence in compar- 
ing a running coupling of the non-supersymmetric model with a parameter of the conformal, 
supersymmetric model, as we discussed above. Secondly, because, even if one could prove that 
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the identification of the A parameter of the SYM theory with the MS 't Hooft coupling is well- 
motivated and should be preferred over other schemes, it is not at all obvious that, at the 
(relatively small) A-values where the YM plasma appears to become quasi-conformal, subleading 
corrections not included in eq. would still be negligible. Finally, because it should be noted 
that the A-dependence predicted by eq. (|3D in a range of couplings comparable to those shown 
in fig. [6] is relatively weak, therefore it may well be that the identification of the SYM A parame- 
ter with the YM running coupling evaluated in some different scheme would yield equally good 
optical agreement; in any case, the systematic and statistical uncertainties involved do not allow 
one to draw any firm conclusion about a "preferred" scheme for this type of comparison. 

Another interesting issue concerns the possibility of non-perturbative T 2 contributions to 
A, which may play a prominent role in the temperature range between the maximum of A/T 4 
(which is attained at 1.1T C for SU(3), and tends to values closer to T c when N is increased) and 
the onset of the regime where perturbation theory is accurate: as it has been pointed out by 
Pisarski [23,24], the SU(3) data in this regime [36] appear to be described well by eq. ([I]). 

Our data seem to confirm that this a general feature of all the gauge groups that we studied in 
this work: fig. [7] shows that the results for A/T 4 (which are plotted with the same normalization 
as in fig. [2]) can be described by: 



In ref. [29], it was argued that this behavior may be related with the existence of a dimension-two 
condensate above the phase transition. 

In particular, the parameters obtained from linear fits according to eq. fjl6f) show a weak 
iV-dependence. However, the precision of our results does not allow us to rule out the possibility 
that the data may actually be described by different, more complicated functional forms (possibly 
involving combinations of logarithms of perturbative origin). 

4.2 Systematic uncertainties, extrapolation to the continuum and to the large- 



Having presented the numerical results from our main ensembles of configurations (which have 
been obtained on lattices with Nt = 5 for the finite-temperature runs), we now turn to a discussion 
of the main systematic errors in our calculations, and to the issues related to the continuum and 
large-iV extrapolations. 

The main systematic uncertainties involved in our calculation include the setting of the phys- 
ical scale, discretization effects, and finite-volume effects. 

To set the physical scale for our simulations, we used high-precision non-perturbative results 
which are available in the literature. In particular, for the SU(3) gauge group, we set the scale 
using the tq parameter, via the interpolating formula given by eq. (2.6) in ref. [133], which holds 
for 5.7 < (3 < 6.92, and thus completely covers the range of points that we simulated. Therefore, 
the systematic uncertainty associated with the scale determination for our SU(3) data receives 
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Entropy density vs. 't Hooft coupling 
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Figure 6: (Color online) The entropy density s (normalized to its limit for the non-interacting 
gluon gas so) versus the renormalized 't Hooft coupling A evaluated in the MS scheme at the 
momentum scale characteristic of the corresponding temperature, as estimated from two-loop 
perturbation theory; note that the data correspond to temperatures which are increasing from 
left to right. The dashed curve is the AdS/CFT prediction for the dependence of s/so on the 
't Hooft coupling, at the leading order in an expansion around the strong-coupling limit, eq. 
intriguingly, if one identifies the 't Hooft coupling of the supersymmetric model with the running 
coupling A of the YM theory in the MS scheme, the numerical results in the temperature regime 
where the deconfined plasma approaches a quasi-conformal regime — see fig. [5] — tend to get very 
close to the curve describing the AdS/CFT model. However, see the text for a discussion of the 
ambiguity in such a comparison. 
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Figure 7: (Color online) Does A receive non-perturbative 0(T 2 ) contributions in the range be- 
tween the deconfinement temperature T c and the onset of the perturbative regime? Do such 
contributions appear only in SU(3), or are they a feature of the large- iV limit too? This diagram 
shows that the data of fig. [5] for the A/T 4 ratio, plotted against (T c /T) 2 , may be compatible 
with a linear behavior. The dashed straight lines are obtained fitting the numerical results to 
eq. ||ISJ in the range (T c /T) 2 < 0.9, with m = 1.64(1), 1.60(2), 1.480(9), 1.67(1), 1.622(9) and 
q = 0.006(3), 0.028(6), 0.037(5), -0.005(3), 0.005(4) for SU(3), SU(4), SU(5), SU(6) and SU(8), 
respectively. 
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contributions from the accuracy limits on that interpolating formula (which are between 0.5% 
and 1%), and from the ambiguity in choosing the dimensionful quantity through which the scale 
is set: different observables would lead to slight, 0(a 2 ), differences in the scale determination. For 
the SU(iV > 3) groups, we set the scale by interpolating the string tension values calculated by 
Lucini, Teper and Wenger [134, 135]; since these results do not extend up to the largest (3- values 
that we simulated, in the latter range we determined the scale following the strategy discussed in 
ref. [138], namely through the three-loop perturbative lattice /^-function [136] in the mean-field 
improved lattice scheme [137], taking lattice corrections into account (note that in ref. [138] it 
was shown that the two methods are, in fact, consistent). The systematic uncertainty associated 
with the scale determination is difficult to quantify, but can be (roughly) estimated by comparing 
the results obtained using different interpolations for log(cra 2 ) as a function of ((3 — /?*) (where /?* 
is a value chosen within the interpolating range). We estimate the systematic uncertainty related 
to the scale determination to be up to 2%. 

To get a reliable extrapolation to the continuum limit, in principle one should repeat the whole 
calculation on a series of lattices of increasingly finer spacing, i.e. at increasingly larger Nt values 
(keeping the Nt/N s ratio fixed, to disentangle from possible finite- volume effects), and perform a 
linear extrapolation in a? to the limit a — > for each gauge group separately. Unfortunately, for 
a set of high-precision calculations like the present one, this task would be computationally very 
demanding. On the other hand, in the literature there are indications suggesting that the results 
from our simulation ensemble A, obtained on lattices with Nt = 5, may already be relatively 
close to the continuum limit. In particular, ref. [96] discusses a recent calculation similar to ours, 
finding that, for SU(4), the continuum limit is reached at Nt = 6. 

To check the impact of discretization effects on our results from the Nt = 5 lattices, we 
also performed a preliminary set of SU(3) simulations on lattices with Nt = 10, N s = 40, that 
we denote as 'ensemble in the following — see tab. [2] for technical information. However, 
the comparison between the results for the trace of the energy-momentum tensor in SU(3), as 
evaluated from the two ensembles, which is shown in fig. El does not give conclusive information 
on discretization effects. By comparing with similar works in the literature [96], we take 3% as 
a rough estimate of the impact of discretization effects (not accounted for by the rescaling of 
the SB limit on the lattice through the factor Rj(Nt)) affecting our main set of results from the 
ensemble A. 

As it concerns finite-volume effects, at the temperatures investigated in the present work one 
generally expects them to be strongly suppressed, due to screening. In ref. [99], we checked this 
explicitly for the setup of the present simulations, and found that further, non-trivial infrared 
corrections to the lattice equation of state [142] are not visible within our data precision at the 
temperatures and volumes considered in the present work. 

To summarize, we estimate that the major sources of systematic uncertainties in our simu- 
lations come from the determination of the scale and from discretization effects. In particular, 
in our plots, the scale uncertainty has an impact on both the temperature axes and on the ther- 
modynamic quantities, which have an explicit scale-dependence through eq. (|13|) . Combining the 
errors associated to scale determination and discretization effects in quadrature, the systematic 
uncertainty in our results turns out to be around 3%. 
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Discretization effects in SU(3) 
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Figure 8: (Color online) Comparison of the results for the trace of the energy-momentum tensor 
in SU(3) obtained from our main set of simulations (ensemble A, corresponding to lattices with 
Nt = 5 and N s = 20, denoted by the small black symbols) and from a set characterized by a 
lattice spacing twice as fine (ensemble Q, with Nt = 10 and N s = 40, denoted by the large violet 
symbols). 
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N s N t 


T = statistics finite-T statistics 


(3 T 
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9000 2800 


6.055 0.804 
6.180 0.978 
6.305 1.170 
6.430 1.383 
6.550 1.613 
6.680 1.903 
6.805 2.239 



Table 2: Parameters of the set of lattice simulations (ensemble O) used to check the impact of 
discretization effects in the main set of results presented in this work. 



Finally, we perform an extrapolation of the results from our ensemble A for A/T 4 to the 
A — ► oo limit. According to the method we used in this work, this observable is evaluated from 
direct measurements of plaquette expectation values, see eq. (fl"3j) . whereas all other quantities 
are obtained indirectly from it, as discussed above. 

In order to do the extrapolation to the large-A limit, we first interpolate our results for the 
trace of the stress-energy tensor for the different gauge groups, according the functional form 
given in eq. (CI) of ref. [42] (redefined according to our normalization conventions): 
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which models the region near the critical temperature using a modified hyperbolic tangent and 
includes T~ 2 and T~ 4 terms to mimic the data behavior at temperatures intermediate between 
the A/T 4 maximum and the onset of the perturbative regime. The fit results are given in tab. [31 
where the uncertainties shown do not include systematic errors. 
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xVd.o.f. 


3 


1.0317(13) 


0.02557(96) 


1.686(14) 


-0.038(24) 


2.30 


4 


1.0154(20) 


0.0082(12) 


1.763(20) 


-0.185(33) 


2.10 


5 


1.0094(98) 


0.0056(10) 


1.729(83) 


-0.205(65) 


1.56 


6 


0.989(37) 


0.009(87) 


1.803(65) 


-0.244(97) 


2.77 


8 


0.9944(46) 


0.0058(51) 


1.720(26) 


-0.165(40) 


6.11 



Table 3: Results of the fits for A according to eq. JT 
uncertainties. 



errors shown do not include systematic 



The parameters thus obtained are then fitted as functions of A 2 including a constant, a 
linear and a quadratic term: 



(18) 
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we note that, for fi, and fa, the 0(N~ 4 ) term can be neglected, therefore for these parameters 
we perform a two-parameter fit, setting qi to zero (and we include the difference between such 
linear interpolation and a quadratic interpolation as a contribution to the overall systematic 
effects). On the contrary, for fa it is necessary to include the 0(N~ 4 ) term in the fit, otherwise 
one would find an unphysical, negative result for the extrapolated value 02- 
The fit results are displayed in tab. d] and in fig. [9j 



parameter 
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Qi 


xVd.o.f. 


fx 


0.9918(20) 


0.361(21) 




0.32 


fa 


0.00895(53) 


-0.217(16) 


3.30(11) 


0.01 


fa 


1.768(36) 


-0.66(41) 




2.16 


fa 


-0.244(46) 


1.73(55) 




1.56 



Table 4: Results of the fits for A according to eq. (fTTj) : for fa, fa and fa, we perform a two- 
parameter fit, including the q and li terms in eq. (117p . setting qi = 0. However, this procedure 
would yield a non-physical, negative result for the fa extrapolation, therefore for this parameter 
we performed a three-parameter fit, fitting Cj, li, and also q^. The errors shown do not include 
systematic uncertainties. 

Finally, from the parameters thus extrapolated, we obtain the curves for the large- N limit of 
A/(iV 2 T 4 ), of p/(iV 2 T 4 ) (by numerical integration of the latter over logT), and of e/(iV 2 T 4 ) and 
s/(N 2 T 5 ) (through linear combinations of the former two quantities). The results are shown in 
fig. [10] (note that the quantities in this plot are not normalized to their SB limits; the latter are 
denoted by the three horizontal bars on the right-hand side of the figure), where we also show the 
large-iV limit of the latent heat L h calculated in ref. [155], namely: L^N'^T' 1 = 0.766(40). 
Our estimate for the same quantity is fully consistent: L^ 4 A^~ 1 / 2 T~ 1 = 0.759(19). In this figure 
we also show our estimate of the errors at six reference temperatures T/T c = 0.9, 1, 1.5, 2, 2.5 
and 3. 

5 Conclusions and outlook 

In this article, we have presented the results of a lattice study of the main equilibrium ther- 
modynamic observables in SU(iV) gauge theories with N = 3, 4, 5, 6 and 8 colors at finite 
temperature, and discussed their extrapolation to the large-iV limit. Our main motivation for 
this study was to investigate whether (and to what quantitative extent) the large-iV limit cap- 
tures the non-trivial thermodynamics of the deconfined plasma in the theory with N = 3 colors. 
This is crucially relevant for AdS/CFT models or holographic QCD models which are expected 
to provide an analytical theoretical description of the strongly-interacting QCD plasma produced 
in relativistic collisions of heavy ions. To address this issue, we have performed a high-precision 
numerical study, based on a fine temperature scan in the phenomenologically most interesting 
region between 0.8T C and 3.4T C . 
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Figure 9: (Color online) Extrapolation of the parameters of eq. (|17p . as obtained from fits to our 
data for the different SU(iV) groups, to the large-iV limit. For each parameter, the result in the 
N — ► oo limit is obtained through a first- (for f%, and f^) or second-order (for $2) polynomial 
extrapolation in N~ 2 , according to eq. (|18p . 
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Extrapolation to the large-iV limit 




Figure 10: (Color online) Final results for the large-iV limits of p/(N 2 T 4 ) (black solid curve), 
A/(iV 2 T 4 ) (red dashed curve), e/(iV 2 T 4 ) (blue dash-dotted curve) and s/(N 2 T 3 ) (green dotted 
curve); the errorbars (including statistical and systematic uncertainties) at six reference points 
corresponding to T/T c = 0.9, 1, 1.5, 2, 2.5 and 3 are also shown. In the normalization used in 
this plot, the Stefan-Boltzmann limits for the pressure, energy density and entropy density are 
respectively denoted by the three horizontal bars displayed on the right-hand side of the plot, 
from bottom to top. The maroon arrow on the left-hand side of the figure denotes the large- N 
limit of the latent heat L^, as calculated in ref. [155], which is consistent with our result. 
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Our results, which agree with previous similar studies [94,96], show that the pressure, the 
trace of the stress-energy tensor, the energy density and the entropy density (per gluon) evaluated 
in all gauge theories that we considered are very close to each other and to those obtained in the 
SU(3) model. Essentially, the main difference between the results we obtained for SU(3) and for 
the other SU(iV > 3) groups is that, as expected on general grounds [156-158], the latter exhibit 
a more strongly first-order deconfining phase transition. 

This led us to compare our results with the improved holographic QCD model recently pro- 
posed in refs. [77-80], which is an AdS/QCD model providing a quantitative description of the 
features of the strongly-coupled regime of the gauge theory in terms of a dual, five-dimensional 
gravity model with a dilaton potential depending on two free parameters. Our results turn out to 
be in very good agreement with the curves from the improved holographic QCD model published 
in ref. [82], as fig. CO fig. El fig. 1 and fig. H show. 

In the improved holographic QCD model (and in other similar models) the running of the 
coupling in the gauge theory is associated with a non-trivial dilaton in the dual gravity theory; in 
particular, the parameters of the dilaton potential can be chosen consistently with the quantitative 
results known for the gauge theory, both in the perturbative and in the non-perturbative regimes. 
This successfully reproduces the non-trivial features observed in the thermodynamic quantities 
considered in this work, including, in particular, their behavior in the temperature range above 
and close to the deconfinement temperature. 

The latter temperature range is indeed particularly interesting, not only from the phenomeno- 
logical and experimental point of view, but also from a theoretical perspective, since it is a 
regime where weak-coupling expansions of the pressure typically show large deviations from lat- 
tice results, and, in particular, fail to reproduce the large values of A observed in numerical 
simulations — see ref. [22] and references therein. In particular, the observation of a possible In- 
dependence of A in the lattice results for SU(3) from ref. [36] led to conjecture the existence of 
contributions which could not be captured by a perturbative approach [23,24]. Our results appear 
to confirm that the T 2 -behavior A is common to all gauge groups that we studied, and thus is 
likely to also persists in the N — > oo limit (see fig. [7|), but we cannot rule out the possibility that 
our data could also be fitted by some more complicated functional form. 

Next, we discussed some issues related to the description of the QCD plasma in terms of 
AdS/CFT models, in particular the large- N limit of the M = 4 supersymmetric Yang-Mills the- 
ory. While the gauge-gravity correspondence [57] provides one with formidable tools to study 
the strongly-coupled regime of this theory analytically, via calculations in the weakly-coupled 
regime of the dual supergravity theory, the supersymmetric theory is not a close approximation 
of real-world QCD; in particular, the former is conformally invariant, whereas in the latter the 
coupling runs with the momentum scale. However, it is possible that QCD may admit an effec- 
tive AdS/CFT description, at least in a regime where it tends to be almost scale- invariant. In 
fig. [5l we have shown that, for temperatures around 3T C , the results for the pressure versus the 
energy density for all gauge groups are already relatively close to the conformal limit e = 3p, 
while the theory is still strongly interacting and far from the Stefan-Boltzmann limit. It is in- 
triguing to observe that, if one identifies the physical 't Hooft running coupling evaluated at such 
temperatures in the MS scheme with the A parameter of the M = 4 model, then the supergrav- 
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ity prediction for the ratio of the entropy density to its limit in the non-interacting gas [122] is 
in very good agreement with the lattice results precisely in the expected temperature range, as 
fig. shows. Although, as we have discussed, such comparison between the Yang-Mills models 
and the supersymmetric theory may be questionable at a fundamental level, prone to scheme- 
dependence, and the observed agreement may be accidental (because the dependence of s/s$ on 
A in the supergravity theory is only known at the first non-trivial order in an expansion around 
the strong-coupling limit), we also note that, if one just compared the lattice results with the 
s/sq value predicted by the supergravity model in the infinite-coupling limit, which is 3/4 [121], 
then in the temperature regime where the thermodynamic observables (in units of the appropri- 
ate powers of the temperature) tend to become almost T-independent the lattice results would 
overshoot the AdS/CFT prediction. This observation may be relevant to identify the parameters 
of AdS/CFT models appropriate to describe the QCD plasma [71]. 

Finally, we have discussed the systematic uncertainties affecting the calculations presented 
in this work (including, in particular, those related to the scale determination, to discretization 
effects, and to finite- volume effects), and performed an extrapolation to the large- .A limit of the 
various thermodynamic quantities; the results — for a parametrization of A according to eq. (j 1 71) — 
are reported in the second column of tab. H] and displayed in fig. [10] (in which we also show the 
estimated errorbars at some reference temperatures). 

To summarize, our results show that the non-trivial features of equilibrium thermodynamic 
observables in the strongly-interacting deconfined plasma are common to all gauge groups in- 
vestigated in this work, and to the large-A limit. Theoretical models based on the analytical 
simplifications and on the dualities that are conjectured to be realized in that limit are thus 
expected to be relevant for the thermodynamic properties of the plasma of SU(3) Yang-Mills 
theory (and, likely, of QCD) at finite temperature. This is particularly important to get ana- 
lytical predictions about quantities (e.g. transport coefficients, observables involving real-time 
dynamics) and/or regions of the QCD phase diagram (e.g. at large densities) where lattice QCD 
calculations are hindered by technical difficulties. 

As for possible future research lines related to this project, it would be interesting to extend 
the present investigation of the thermodynamics of hot SU(-A) gauge theories by looking at 
other observables, for instance screening masses, and comparing the results with predictions 
from the improved holographic QCD model [77-80] or from AdS/CFT [160]. The investigation 
of observables related to non-equilibrium thermodynamics in the large-A limit would also be 
interesting, but is expected to be much more demanding from the computational point of view. 
On the contrary, another interesting research direction would be to perform similar studies for 
models defined in lower dimensions [168], and/or based on different gauge groups, not necessarily 
relevant for the large- A limit, for which powerful numerical algorithms are available [169-176]. 
This could lead to a better identification and understanding of the universal and of the model- or 
dimension-dependent features of strongly interacting systems at finite temperature, and possibly 
shed new light onto the fascinating mysteries of the QCD plasma. 
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